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^-h Abstract 

p^j In the context of complex WKB analysis, we discuss a one-dimensional Schrodinger equation 



h z d z x f{x, h) + \Q{x) + hQ x {x, h)]f(x, h) = 0, h-+0, 



where Q(x), Q\{x,h) are analytic near the origin x — 0, Q(0) = 0, and Q\(x,K) is a factorially 
divergent power series in h. We show that there is a change of independent variable y — y(x, h), an- 
alytic near x = and factorially divergent with respect to h, that transforms the above Schrodinger 
equation to a canonical form. The proof goes by reduction to a mildly nonlinear equation on y(x, h) 
and by solving it using an appropriately modified Newton's method of tangents. 

1 Introduction 

\q Context. The stationary Schrodinger equations for an anharmonic oscillator 

^ -h 2 d 2 J(x,h) + Q(x)f(x,h) = 0, (1) 

o 

where h — > 0+, f(x, h) is an unknown function, and V(x) is, say, a polynomial with 

Q(0) = Q'(0) = + Q"(0), 

^ is one of the basic problems in Mathematical Physics. It is easy to find formal WKB solutions of (fTl) 

H by substituting an ansatz 

f(x, h) ~ e s ^/ h (a {x) + ai(x)h + a 2 {x)h 2 + ...) (2) 

and recursively solving for new unknown functions S(x), ao(x),ai(x), .., yet analytic properties of this 
expansion are quite subtle. Even in simplest examples, the terms of the series ao(x) + a\{x)h + 
a>2(x) h 2 + ... have singularity at the origin and diverge factorially (i.e., are Gevrey series, sec. J3|) away 
from zeros of Q(x); a lot of effort is being put into proving Borel summability divergent series of this 
origin. 

The singular behavior at the origin of the coefficients dj(x) should not be all that uncontrollable: 
the closer x gets to 0, the more the equation (fTl) resembles the Schrodinger equation of a harmonic 
oscillator whose eigenfunctions can be expressed in terms of well-known cylindro-parabolic functions. 



And in fact, in [AKT09J it is shown that one can find a change of coordinates y = y(x,h), where y 
is a no more than factorially divergent power series in h with coefficients holomorphic functions of x in 
a full neighborhood of the origin that reduces (fij) to a harmonic oscillator. The proof goes by reducing 



the problem to a mildly nonlinear equation (case M = 2 and Qi(x, h) = of our (13), see also (14)), 
recursively solving it for Eo tn ,T n (z) where Eq = ^2 n>0 Eo !n h n , T = ^2 n>0 T n (z)h n and analyzing the 
growth of |-Eb,n|> |^n(z)|- Much more is contained in the rich and beautiful paper [AKT09] . but here 
we will concentrate just on this aspect. 

Motivation. Our interest in the Witten Laplacian [Gllj led us to study the paper [AKT09] 
with the task to extend its results to equations of slightly more general type than (fiT), namely: 

- h 2 d 2 J{x, h) + (Q{x) + hQi(x, h))f(x, h) = 0, (3) 

where Qi(x, h) = J2 n >o Qi,n(x)h n is a factorially divergent series in h, Q{0) = Q'(0) = ^ Q"(0); in 
particular, to construct a coordinate change y = y(x,h) in a full neighborhood of x = that would 
deform (T3J) to a harmonic oscillator. By the time we asked all why and why not questions about the 
original proof of [AKT09J , we ended up with their proof repackaged in the form of Newton's method 
of tangents in the spaces of factorially divergent series. This version gives us the desired generalization 
automatically. Just as automatically we have obtained a generalization to an arbitrary order of zero 
ofQ(x). 

The main result is: 

Theorem 1.1 Suppose that in the equation (pi). - 

a) all functions Q(x), Qi n (x), n > 0, are analytic on a common neighborhood U C C of x = 0; 

b) Q(0) = Q'(0) = ... = Q( M - l \0) = + Q( M )(0), M > 1, 

c) there are constants r, Cq such that sup^g^ |Qi, n (a ; )| < C§T n n\. 
Then there exist: 

1) series y = y(x,h) = J2n>oyn( x )h n , where y n {x) are holomorphic functions on a common neigh- 
borhood [/'cC of x = 0, 

2) series Ej{h) = £n>o E j>n h n , E j;Tl G C, for 0<j<M-2, 

3) series ifi = ip(x, h) = X^n>o^™( x )^™' ^nix) analytic on U' ; 

4) constants C' Q , r' 

so that 

sup \y n {x)\ < C' T /n n\, sup \ipn(x)\ < C' T /n n\, \E j>n \ < C' T /n n\, 
xeU' x&J 1 

and the change of independent variable y = y(x, h) followed by rescaling of the unknown function 
f(x,h) = ip(x,h)g(y(x),h) in ([3]) transforms this equation to 

-h 2 d 2 y g{y, h) + (hE (h) + hEh(h)y + ... + hE M . 2 (h)y M - 2 - y —)g{y, h) = 0. 



Plan of the paper. Section [2]is devoted to formulating our problem in terms of a mildly 
nonlinear equation ( |13[ ). Sections p] and ^introduce definitions of functional spaces used in the rest of 
the paper. The section [5] discusses inversion of the linear operator which is the dominant member of 
equation (13). Section p\ contains no precise definitions or results but explains the intuition based on 



the Newton's method behind our argument. Section [7] contains quite general results about analytic 
operations with factorially divergent expansions. The main result is proven in section |8] 

The contribution of this paper is, in our understanding, as follows. Theorem |1.1| is a new 
generalization of |AKT09l Theorem A.l]. In the vast literature on Gevrey solutions of nonlinear 
differential equations we have not seen any applications of the Newton's method of tangents, so this 
way of approaching the problem is probably also new. The intermediate result of section [7] should be 
well known, but at the moment we do not have a specific reference. 

2 Formal reduction to the normal form. 



Following the ideas of |AKT09l §1 and Appendix A] whose treatment corresponds to ours in case 
M = 2, Qi = 0, we consider an equation 



h 



2 d 2 f(x,h) 
dx 2 



+ (Q(x) + hQi(x,h))f(x,h) = 0, 



(4) 



where Q(x) is an analytic function in a neighborhood U C C of the origin, Q(x) = cx M + 0(x M+1 ) 
as x — > for some c / 0, and Qi(x, h) = Yl n >o h n Qi,n(x) is a formal (for now) power series in h with 
coefficients Qi t n(x) analytic functions on a common neighborhood of x = 0. 

We seek to simplify Q by an invertible change of independent variable y = y(x, h) = ^2 yj(x)h? , 
which is (for now) a formal power series in h with coefficients yj(x) holomorphic functions in a common 
neighborhood of 0. Then Q becomes 

lr I -4 -— 7 i — i —I + ( — ) (Q(x) + hQi(x, h))f(x) = 0. (5) 



, i d 2 f d 2 x / dx\ l df 
\dy 2 dy 2 \dyj dy J \dy J 



Let us further replace the unknown function f(x, h) with a product ij){x, h)g(x, h), where g(x, h) 
will be the new unknown function and ip will be conveniently chosen: 



-h 2 - 



dx 



dy dy 



dy 2 dy dy dy 2 



dy 2 \dy 
The following choice eliminates the summands containing jS-: 



(Q(x) + hQi(x, h))i/)(x)g(z) = 0. 



dy- 



^' fe)=eXP U/^O ^ 



The equation Q is thus reduced to 



^d 2 g 

dy 2 




dx Y(Q(x) + h Ql (x,h)) + h 2d2x [dxVl { ^ ,3y) ^(^/^ 



dy 2 \dy J 



i> 



h 2 



p 



(6) 



9 = 0- (7) 



We will now achieve that the expression in braces in ffh takes a form ^7=0 hEj(h)y 
power series Eq(K), ...,EM-2(h). Let us spell out this condition: 



yM 

J — ^j- for 



(()(,')-!-/,(>!. { .r.],))+h 2 [^ 



d 2 x ( dx\ l (dil)/dy) (d 2 ip/dy 2 ) 



dy 2 \dy J 



W 



V 



dy 
dx 



2 Af-2 

3=0 



,o_v_ 



A one-page elementary calculation allows us to rewrite the /i 2 -term on the left and obtain an 
equation 



j \ 2 M-2 M , 2 

Qi.r) + h(hU.h) = (2) (^2 hE^hW - y —) - - 



dx I 



3=0 



d^y f dyy 1 _ 3 / tfy\ 2 ( dy 



dx 3 \dx J 2 \dx 2 J \dx 



(8) 



Following |AKT09J . in the equation (J8J) we replace the independent variable x by z, where 

2 
M+2 



M 



z = A(Jq y/-Q(t)dt) «+ 2 . It will be convenient to choose A = (M + 2) 

Quite generally, change of independent variable can be performed using formulas 

dy dz dy 
dx dx dz ' 



d 2 y _ dyd 2 x f dx\ 3 d 2 y / dx\ 2 
dx 2 dz dz 2 \dz J dz 2 \dz J 



d 3 y d 2 yd 2 x ( dx\ dyd 3 x ( dx\ n dy ( d 2 x^ 

dx 3 dz 2 dz 2 \dz J 



dz dz 3 \dz J dz V dz 2 



dx 
dz 



+ 



d 3 y f dx 



dz 3 V dz 



The equation Q becomes: 

Q(x) + hQ 1 {x,h) 



, \ 2 / , \ 2 M-2 



dz 



dx J 



3=0 



M 

4~' 



2 



d 2 yd 2 x f dx\ dyd 3 x ( dx\ dy fd 2 x 

dz 2 dz 2 \dz J dz dz 3 \dz J dz \ dz 2 



)(£r+SK£)> 



dy d 2 x f dx 
dz dz 2 \ dz 



+ 



d 2 y f dx 



dz 2 \ dz 



dz 



dy 
dz 



(9) 



Multiplying both sides by (^) and using 



we get 



dz i 



dx\ 

~d~z) 



-A 



-M-2 



M + 2 2 z M 
2 } Q{x) 



-A 



2 r' +h (di) Q ^ h) 



-M-2t M + 2 ^2M ■ , (te\ 



2 



-2 



dv\ 2 M ~ 2 i, M 



dz 



3=0 



d 3 xdz ^ 3 fd 2 x y (dz_\ 2 ^ <$y_ (dyY 1 _ 3 ( f]L\ ( ^ ' 

dz 



dz 3 dx 2 V dz 2 J \dx J dz 3 \ dz 



To make the structure even more transparent, put 

-2 



Qi{z,h):=Y,h n Qi t : 



n>0 



w =- (I 



_(z(»,/i) -/i 



2 \dz 2 



d?xdz 3 (d 2 x^ 2 (di 



dz 3 dx 2 V dz 2 J \dx 



(10) 



Then ([To]) becomes (recall that A = (M + 2) m+2 

v M K l 



Z i + HQ l{zM = (§ 



£ E iW 



y 



3=0 



dz 3 \dz 2 l dz 2 ) [dz 



• (11) 



From (jl Xp we can see directly that the h° term in y(z, h) equals to z. Thus, we can replace the 



unknown y of (11) by T as follows: 



y(z, h) = z + hT(z, h). 



In terms of T(z, h) the equation (11) becomes (with primes denoting 4-] 



M 



M-2 



+ hQ 1 {z,h) = {l + hT'yC£^E 3 {h){z + hTy 

3=0 

which after cancellation of h° terms becomes 



(z + hT) M , h 



hT'" 3 (hT 



n \2 



1 + hT 2(1 + hT') 2 

(12) 



M-2 



y M _ ( M \h^- 1 z M -^Ti j - 

F(T,E ,...,E M „ 2 ) := J2 E 3<h)^ + hTy - ^~ lK ^ 



+ 



j=0 



-T (2 + hT') 



M ~ 2 \ U2 rplll 0,3 rpl,2 

(13) 



Proposition 2.1 The relation 

F(T,Eo,...,E M -2) = 

seen as an equation on formal power series in h Ej(h) = Yln>o h n Ej,n, Ej tU G C, and T(x,h) = 
Y2 n >oT n (z)h n , has a solution with T n (z) holomorphic in any connected neighborhood of z = on 
which all Qij(z) are holomorphic. 



Proof. A simple recursion with respect to the power of h can be set up similarly to [AKT09, 
Th.1.1, Rmk 1.1] once we notice that 

F(T, E , .., E M -2) = E + ... + E M -2Z M ~ 2 - Izf(zfT)' - Q^z, h) + 0{h). (14) 



□ 



in h. 



It is the goal of the rest of the paper to estimate the growth of term of the obtained power series 



3 Norms on Gevrey series 

The content of this section [3] is classical. 

Let U be an open subset of C, let P(z, h) = Ylk>o h k Pk{z) be a formal power series in h with 
holomorphic coefficients Pk(z) G 0(U). 

We say that P(z, h) is Gevrey-type on U if for any compact subset K C U, there are constants 
Mki Pk such that 

sup \p k {z)\ <M K p k K k\. 



Consider the space of those P(z, h) = P(h) which do not depend on z, endowed with the norms 
which depend on a parameter t > 0: 

No(P,t) = Y,^t n - (15) 

fc>0 

An easy calculation establishes: 
Proposition 3.1 Suppose P(h),Q(h), and t are such that No(P,t), No(Q,t) < oo. Then: 

N (PQ,t) < N (P,t)-N (P,t), 
N (P + Q,t) < N (P,t) + N (Q,t). 

Let V* be the subspace of those P(h) for which No(P,t) is finite. Analogously to the proof of 
completeness of £ 1 (C), one verifies that V* is a Banach space. 

Following |Sch| p. 15], for a Banach space A and a number p > denote by A( / o) the space of all 
formal series {g = J2j>o a i T "'! a j e A} such that 

Nip := ElMI^' ( 16 ) 

i>o 

is finite. The space A(/3) endowed with || • \\ p is a Banach space (which can be proven analogously to 
the proof of completeness of the ^ 1 (C) space). If A is a Banach algebra with ||a6|| < ||a|| ||6|| (e.g., 
if A = V*), then given another such series g = Sj>o%" r i' we bave ||gp||p < \\g\\p \\g\\p for the usual 
product of power series. 

Clearly an element of V*(p) for p, t > gives rise to a Gevrey series P(z, h) for z in a small disc 
around the origin. 

The following Lemma replaces the Cauchy integral formula when it comes to estimating the 
norm of a derivative of ^ • ajT 3 : 

Lemma 3.2 If g = ^2 n>Q a n T n G A(p) and we let g( k > := ^2 n>0 a n n(n — l)...(n — k + l)T n ~ k fork > 1, 
then for any e, < e < p, we have 

II (k) n ^- ""• n n 

\\9 'llA(p-e) < ^ll5llA(p)- 



Proof. Indeed, by the binomial formula 

-k 

k 



J £ \K\Hn - l)...(n -k + l)(p- e)*"* = £ \\a n \\ (?) (p - e) n - k e k < £ 

' n>0 n>0 ^ ' n>0 

□ 



4 Functional analytic setup 



Let us specify the functional spaces between which the nonlinear functional F given in (13) will define 
a continuous map. 

Notation 4.1 Let us fix for the rest of the paper to > 0, po > in such a way that Qi(z, h) defines 
an element of V T °(po)- 



If necessary, shrink po to be < 1/M, this will be used on p 21 



We define the following spaces for 0<s<1,0<£<to: 
X\ := {(E ,...,E M „ 2 ,T)E Vt x • • • x V* xV ^ o( ^ + s) ) : T,T',T",T"' € V*( ^-^ )} (17) 

M — 1 factors 

endowed with the norm 

M-2 3 jj 

11(^0,-", -^A/-2,r)|| S)t := 2^ \\Ej\\yt + 2 , || J7j T Kt( Po( - 1+a) y ( 18 ^ 

i=o i=o J 

(the reason we do not want to consider arbitrarily small neighborhoods of z = is the appearance of 



the radius of a neighborhood in the denominator of the estimates in Lemma 5.1); 



y « :=vl( Mi^) ); m 

z> := {(£„,..., £ M _ 2 ,r) e y' x ■■■ x V, xV(g^ti)) : r,r'6V'(«fl_t£> )} (20 ) 

M - 1 factors 

endowed with the norm 

M-2 1 „■ 

||(£?0, ...,E M -2,T)\\s t t ■■= 2^ ll-^jllv* +2^ \\j-j T \\ v t(PQil±£ly ( 21 ) 

j=0 j=0 

Clearly, for every s,t such that < s < 1, < t < To, F defines a continuous and even 
analytic map F* from a neighborhood ZY* of the origin in X\ to the space y\; for concreteness, Wj = 
{(Eo, ...,Em-i,T : t||T"|| v t( s ) < 1}. 

The spaces iT* are auxiliary and will be used in the proof later on. 

For any fixed t, we have the following properties, of which 1) and 2) are trivialities and 3) follows 
from Lemma 13.21 



1) For s' < s, there are inclusions of norm < 1: X l s ^-> X*,, y\ <^ }>*,, Z>\ *-> <2y- This is 
expressed by saying that Xf •., y\\, Z\ \ are scales of Banach spaces 

2) The obvious map X l s — \ Z\ has norm < 1; 

3) For every s' < s, there is a map Z l s — > X l s , so that 



Xl\\ < 



B 



(s - s>) 



'\2' 



(22) 



where B depends on po but not on t,s,s'. 



5 The dominant term of the equation F(T, Eq, ..., Em -2) = 



In this section[5]we are studying the h° component of the equation F(T,Eq,...,Em-2) = 0, cf. (13). 
The following lemma is inspired by |AKT09( Lemma A. 3]. 

Lemma 5.1 Let M > 1 be an integer, v(z) a holomorphic function on A = {z : \z\ < r^} with values 
in a Banach space B, and consider the following equation for an unknown holomorphic function 
u(z) : A — > B and constants Eq, ...,Em-2 G B: 



z M d M M , 
2 dz 4 



E + zEi + ... + z M ~ z E M -2 + v(z). 



(23) 



Then (23) has a unique solution (u(z),Eq, ...,Em-2) and moreover for any r, < r < r$: 



\Ej\\b < ^IMIfl(r), j = 0,...,M-2; 



\ U \\B(r) < ZMZi\\v\\B(ry, 



du 



dz 



< 



B(r) 



^\\ V \\B(r)- 



(24) 

(25) 
(26) 



Proof. Rewrite the equation (23) as follows: 



1 



z^(z^u(z))' = E + zE x + ... + z M ~ 2 E M -2 + v(z). 



(27) 



If v(z) = vo + v\z + V2Z + ..., then (27) can be rewritten as 

Af-2 



1 Af-2 oo 



■_M 



i=o 



j=M-l 



(28) 



If u(z) is to be holomorphic, the RHS of (28) should not have any terms with z- , hence 



-Ej, ]fj-M.<-l. 



Inserting this, we have 



oo 



l -(zfu(z))'= Y, zi-f(E j+Vj )+ JT * j ~ fv i> ( 29 ) 

jeNo; ¥—l<j<M-2 j=M-l 



or 



■^ „-_ji/f_i_i /„ n ^ ^ £ 



iSNo; f -l<i<M-2 J 2 j=M-\J 2 T 

We conclude that the solution is necessarily of the form 

Ej = -Vj, j = 0,...,M-2, (31) 

oo 

«(*)= E Hbr^~ M+1 - 

j=M-l J 2 "+" X 

As j - 4f + 1 > 5 for j > M - 1, we have 

oo oo 

K*)|| B (r)<4 E \H\bt j - m+1 <^- m Y hiW^^^-^Hz)^). 

j=M-l j=M-l 



Further, 

j=M-l J 2 "+" 1 

hence, as the fraction is always < 1, 

[|u'(*)b (p) < 2 Y hj\\BT j - M < 2r- M \\v(z)\\ B(r) . 

j=M 



From (31) we obviously get (24) . □ 



6 Newton's method 

This section [6] contains no precise results, it will not be referred to except for purposes of intuition. 

6.1 Newton's method in Banach spaces. 

Newton's method of tangents for solving scalar nonlinear equations appears in almost all elementary 



calculus textbooks. In this section 6.1 we review the Newton's method in Banach spaces follow- 
ing [KaAkl Ch. XVIII]. 

Suppose X, y are two Banach spaces, x^ G X, U C X is an open neighborhood, T :U — > y a 
continuous map admitting two continuous and bounded Frechet derivatives. Suppose that for all x 6 U 
the inverse of the Frechet derivative, (dF x )~ 1 : y — > X is uniformly bounded, i.e. || (rfi 7 ^) - x || < A. 

9 



By Newton's method we mean the following iterative procedure of finding a zero of T: 

1) x(°) is given; y(°) = F(x^); 

2) for k > 0, define w^ = -(dF a . (fc ))- 1 (y( fc )); 

3) x( fc+1 ) = x^ fc ) + uA fc ); assume or prove that x( fc+1 ) £ U; y( fc+1 ) = F(x^ k+1 ^). 

Assume that x* is such that F(x*) = 0. Then 

y (k+i) = F ( x (fc+i)) = F(x (fc) ) + [dF (xfc) ](u; (fc) ) +0(||u; (fc) || 2 ) = 0(||y (fc) || 2 ) 

S v ' 

=0 

which is smaller than ||j/ '|| if y' > was small already. So the method converges if x^> was sufficiently 
close to x* and, loosely, far enough from the boundary of U. Formalization of these conditions can be 
found in [KaAkl Ch.XVIII, §1.5]. 

6.2 Newton's method leads to Gevrey expansions 



In this subsection 6.2 is to motivate that factorially divergent or, for brevity, Gevrey expansions natu- 
rally arise when we attempt to apply Newton's method for solving our equation F(Eq, ...,Em-2i T) = 
0. 

The Newton's method requires inverting the Frechet derivative of F. Let us for simplicity 
consider the case M = 2 and try to invert the Frechet derivative of F at the point (Eq,T) = 0. If 
(£o,T) is a tangent vector, then 



<lir\A!:\,T)-f ti -^T'-p--^h 2 r". 



z —' z — - 1 • 2<r"' 



Then dFo as an operator X s — > y\ can be written as G + h H, where H is bounded and 

g ■. (£ ,T) ^ £ - Z ^r - ^t ■. xt^yl 



is inverted by an operator K : y s — > X s i, s' < s, of the norm < , _ c ,-, 2 , cf. lemma 5.1 and (22). Let 
us attempt to write the inverse of dFo as 

(dFo)" 1 = ^2h 2n (-l) n (KH) n K 

n>0 

where {KH) n K is 

y\ A x\ x 4 y si 4 xt 2 4 ... 4 3i 4 y s , (32) 

and where we can think of s n < ... < si as arbitrary numbers subject to s' < s n and si < s. Choosing 
si, -,s n to maximize the product for (s — si)(si — S2)---(sn — s') for fixed s, s', we can prove an estimate 

\\{KHYK : ^N^|< /^r 

for some new constants c', c", but no dramatically better estimates are available. Thus, given v € 3q, 

we can hope to represent (dFo)~ 1 v £ X* by the following expansion whose convergence needs to be 

discussed separately: 

r>\C n 
(dFo)- l v = Y, h n u n , s , \\u n J x t < C'j^—-, (33) 

n>0 ^ ' 



with constants C , C independent of s. This shows that expansions of type (33) appear naturally in 



our problem and motivates our decision to formulate intermediate results in terms of them. 

10 



7 Calculus of Gevrey expansions 

7.1 Definitions 

By a Gevrey expansion we mean, vaguely, a formal expansion ^ n h. n u n , where the straight-font h is a 
formal variable, and u n are elements of some Banach space such that ||tt|| < coc n n! for some constants 



Co and c n , compare (33). We choose however to avoid this notion in the mathematically precise body 
of the article. Instead, the following seems to be a more convenient terminology. 

In this section [7| let X s , y s , < s < 1 be two arbitrary scales of Banach spaces, i.e. we suppose 
that for s' < s there is an inclusion X s <—i X s i, y s <—i y s i of norm < 1. In this section [7J X s ,y s have a 
priori nothing to do with X*, y g in the rest of the paper. 

Suppose that for every s there is a map F s : X s D U s — > y s . We will say that the family F = (F s ) 
is compatible with inclusions if: 

a) for s' < s, the map X s '—¥ X s i maps U s into U s r, 

b) for s' < s, {y s ^ y s >) o F s = F s > o (X s -> X s ,). 

For a family of linear maps F s we always think of IA S as equal to X s . 

Analogously we define what it means for a family of maps (G SjS / : X s — > y s ') s '<s to be compatible 
with inclusions. 

We denote by X\ the inductive limit of X s for s < 1 : 

X X '■= {(g s )o<s<i ■ 9s G X s , and for all s < s' (X s > ^ X s ){g s >) = g s }. 

We assume that on every X s there is an action of a linear operator h : X s — >• X s , compatible 
with inclusions, satisfying \\h n : X s —> X s \\ < — u and similar for y s , for some fixed number r > 
independent of s. We will abuse notation by writing \h\ instead of r. 

We call a linear operator G s : X 8 — )■ y s (or a family of such operators) h-linear \tG s oh = hoG s . 

We finish this subsection by stating the following combinatorial inequality: 

Lemma 7.1 |AKT09[ Lemma A. 4] The following inequality holds for all positive integers j, k satis- 
fying k < j: 

Y, ii!j2!...J*!<4*- 1 (j-* ! + l)!. 

ji + ... +jk = j 
ji,...,jk > i 

7.2 Function evaluated on a Gevrey expansion gives a Gevrey expansion 

For k £ Z>o, we call fifc : X s — > y s an h-k-lme&r map of norm < 1 if flk( v ) = ^k(v,v, ...,f) where 
Ofc : X s x X s x • • • x X s — > y s is symmetric, satisfies ||n(vx, ■■■,Vk)\\y s < ||wi||^ s H^Ha's ••• 11^11^ an d, 
moreover, Q k (...,Vj- 1 ,hvj,Vj+i, ...) = hQk(---,Vj-i,Vj,Vj + i, ...). 

Lemma 7.2 Let f : X s — >■ y s be a map compatible with inclusions, Xq € X\; a € M>o- Suppose that 
for any w G X s such that \\w\\ s < — 

y s 3 f(x + w) = f(x ) + Y,a j ^ 3 (w), (34) 
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where Qj jS : X s — Y y s is an h-j-linear map of norm < 1, compatible with inclusions. 
Let for all n> 1, g n = (g r ^ s ) G X\ are such that 

C n n\ 



\9n,s\\x s < Co 



Then there are H n = (H n>s ) G jDi, n > 1 satisfying 



(1-*)"' 



|i?n, S |k</3(C «)- 



C n n! 
1 - sV 



/or /3(t) = te , suc/i f/iai i/ 



t/ien the series 



converges in X s and 



J— 1 — < 1 and > - 1 — — < — — . 

i-8 4i (1_s)n Coa ' 



n>l 



n>l 



f(x + W s )-f(x )=J^h n H n , a . 

PROOF. Motivated by a "formula" 

" f(x + W s ) - /(s ) = ]T a fe fi fc I ]T fc ni g ni , s , ..., ]T fc n * 



m fc ,s 



fe>l 



,n,>l 



n k >l 



fe>l m>l 



/ j »'fe (,5ni,S) ■■•) 9nk,s) 



ni + ... + rik — m 



n 1 ,...,n k > 1 
to which we will give an analytic meaning in a moment, consider 



flm,s ■ — / j 



(Ik 



/ _, "'k \9ni,st ••■■> 9n k ,s) ■ 



fc>l 



m + ... + n k —m 
m,...,n k > 1 



We have, with help of lemma |7T 



\H m , s \\y s <Y,(C a) k - 



c- 



k>l 



(i- s y 



>, nil...n k l < 



iii + ... + n k = m 
m, ...,n k > 1 



(35) 



(36) 

(37) 
(38) 



(39) 



± 713^ ECO^^m - k + 1)! < (Coa) £(0,a)*4^ ^ ^"^V 

v ; k>\ ' v ; fc>0 ' v 



C m m\ 



compare [AKT09L (A.50)]. 

Then, assuming p6|) and hence 



h n C n 



{\- s y- 



— hi Yl m >ih m H m ,s as well as (37) are absolutely 



convergent, so /(xo + W s ) can be written in terms of Taylor series (34), and hence (38) holds by the 
algebraic manipulations of (17. 2|). □ 
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7.3 Increment of a function is well approximated by the first derivative 

Let us now give a meaning to the formula f(x + h b v ) = f(x) + [f'(x)](h b v) + h 2n O(v 2 ) for x = xo + g 
in the context of Gevrey expansions. 



Let Mt) = E k >o {J ^~t^\ 2 (t) = £ fe>1 *£*. 



Lemma 7.3 Suppose f : X s — > y s , xq, a,Qj are as in (34), and g n E X\, n > 1, satisfy (35). Assume 

C>1. 

Then for any integer b > and any sequence v k E X\ with \\vk, s \\x s < -A ? -,_ \k > there are : 

a) a sequence of elements Wj E j^i with \\wj, s \\ < gr-ffi, {Cqo) £'_ ^ satisfying the property: 
under condition 

\h\C ^ h n C n 

< 1 x 



^- h n C n 1 ,„ . 

L £<r^<cw WO) 



f 1 -"' ' ^i* 1 - 8 )" c » a 

we have two vector defined by an absolutely convergent sums 

oo 

v s ■= h b Y^ h k v kjS E X s , g s := ^ h k g k , s E X s 
fc=i fc>i 

cm(i suc/i £/ia£ 

f'(xo + g s )(v s ) = ^2ja J nj(g s ,...,g s ,v s ) = ^hPwj, s in y s ; (41) 

6J a sequence of elements Uj E J^i, j > 1, wit/i 



«j>b B </32(a(C + ^)) 



(1 - s)i 



satisfying the following property: 
under conditions 



oo 



\h\C ^ /i n C 

< 1 



*' ^(l-s)n < (Cn+A)a (42) 



we have 

f(xo + g s + Vs)-f(xo + g s )-[f'(x + g s )](v s ) = h^^Vu^ iny s . (43) 



Proof. 



Preliminary remark. If f is a formal in h written as v = h Xm>i h v k as in the statement of 
the lemma, then it can also be written v = J2 m >b+i h m v m with v m = v m _& and hence satisfying 

|| u m,s|| ^ ^ n_ s )m-i) — -^(l-s)™ II O ^ 1. 

Proof of a) . Rewrite ( |41[ ) as the following "formula" whose analytic meaning will be clarified in 
a moment: 

« £ Y. J^h k ^- +k ^n j (g kl ,...,g kj _ 1 ,v k .) = £>"u; n ." (44) 

J>1 fei,.. ) fcj_i,fcj , >l n>l 
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Thus we define 



Wn = E E 3 C o AaJ ^j(9k 1 ,-,9k^ 1 ,Vk j 

i>! ki,..,kj-i,kj >i 

fci + ... + fcj-i + fcj = n 



and with help of Lemma 7.1 



|Wn,.|k < E E ^ C ^ " ' 

fci + ... + kj = n 



(1 - «)' 



<E^^~^^7i 



L.C-^Cn-^l)!..^ J E(/+1)4 / M /,' I ™ 



i>i 



(1-a)' 



< 



Co lf">o 



J'! (1 - *) r 



Under condition (40), the equation (44) is literally true in y s , both sides being absolutely 



convergent series in which we can change the order of summation, hence (41) 
Proof of b). Let us now formally write 

f(x + g + v)- f(x + g)- [f'(x + g)](v) " = " 



(use symmetry and polylinearity of Clj 

oo 

= E« J E hkl " 

j=2 ki,...,kj>l 
(since Vk = for k = 0, ..., b ) 



E aJ l n j(9 + v ) -^j(g) -j&j(g,-,g,v) 



.0=2 



E I J %(£fcl>-"> 5fc<r,% CT +l,-",% 



E« J E 

J= 2 ki,...,kj- 2 >l; 
kj-i, kj > 6+ 1 



/jfci+-+fcj 



.(7=2 



Ei^j^^i'- 



5 S'fco- 5 ^fccr+1 > •■•' ^fc] 



OO 

E<^ E 

J= 2 ki,...,kj-2>l; 

kj-i, kj > b + 1 



l l k 1 +...+k J 



.(7 = 2 



E ( )^i(5fci>— >3fe -)«fc .+i>->«fcj_2> ,y fej-i-6> v fci-6) 



= E a " E h k ^- +k ^ +2b 

j=2 k u ...,kj>l 

This motivates the definition: 

00 

u « = E a ' E 

i= 2 ki,...,k 3 ■> 1 

fci + ... + fcj = n 



.(7=2 



E ( )^i(fffci'"-'5 , fc [T ,Wfc a+1 ,"-,%_ 2 ,WA :j _ 1 ,UA :j 



.(7=2 



E ( )%(fl r fcl'-'fl , Jfco->^ -+l>-'^-2) W fej-l> W Jfcj 
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We have: 



u 



oo 
n,s\\y s <^2<* 3 

k\ + ... + kj — n 



E 



E 

.(7 = 1 



a) (l-s) n 



ki\...kj\ 



< 



00 fin [ °° I 

< ^^(Co + ^^-Hn-i + l)!-— - < ^^(Co + ^4^ 1 --— 



C n n! 

:i-sK 



We conclude that (43) holds under condition (42) in the same way as we did in part a). □ 



7.4 Inverse of a regularly perturbed linear operator 



In this subsection 7.4 we will generalize the familiar statement that if an operator A has a bounded 
inverse, then its small perturbation A + B has a bounded inverse A _1 (l + BA _1 + BA X BA x + ... 
if ||B||||A _1 || < 1. We show that a similar statement holds if B is replaced by a factorially divergent 
series in h with operator coefficients. 



Lemma 7.4 Suppose given a sequence of h-linear compatible with inclusions operators: 

n\C n 



7~L ns '■ %s ~ * 34 : %n,s — ^i" 



l-s) n ' 



n > 1, < s < 1. 



(45) 



Suppose that we also have h-linear compatible with inclusions operators J- s : X s — > y s , with inverse 
Qs '■ 34 —> X s , \\Q S \\ < A, which is also h-linear and compatible with inclusions. 
Then there exist operators, h-linear and compatible with inclusions, 



C n n\ 



Cn,s :y s ^Xs : \\C n ,s :y s ^X s \\< max{A, (C^V^^W -, n > 0, < s < 1 

(1 — s) n 



(1-8)' 



satisfying the following property: 
with notation "Hq := T , 



Y^ U j>s C m , s - J .^ . f 



0, i/JV>l; 



mj>0; j+m=N 

PROOF. Motivated by the "formulas" 

" (j- + n)- 1 = ((1 + uf- 1 ) o F)- 1 = g o (1 + my 1 = Y,(-iyg(ngy 



(46) 



(47) 



fc>0 



and, for k > 1 



«(ng) k = J2h m Y. Hm, s gn m , s g...n nk , s g n , 



m>l 



rij + ... + n k = m 
rii, ...,n k > 1 



whose analytic meaning will be clarified shortly, we put 



i^k,m,s •" 



E 



W ni ,s5^n 2 ,s5-^Jit5. 



ni + ... + n k = m 
m, ...,n fe > 1 
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We have 



He 



k,m,. 



< 



£ 



m + ... + rik = m 
rn,...,n k > 1 



< (^f- 



C" 



(l-«)' 



(1 - s) n i (1 - s)"2 '" (1 - s )n k - 



^2 n 1 \...n k \ < 



ni + ... + iik — m 
m, ...,n k > 1 



(use lemma 7.1) 



Put for m > 1: 



f~rm 

<{C x Af^—^A k ~\m-k + 1)!. 

J~-m,s : = / { J-J yi^k,m,s\ 



fc>l 



then 



c 



\C m ,s\\ < (CiA) J2(C 1 A) k ——4 k - 1 (m -* + !)!< 



fe>l 



C 



C m m! 



< (^^) 2 7i^ £(^4*v-*')! < (CU) 2 ^^* 



A) 



fc'>0 



(1-a)' 



Put further Cq s = Q. 



(48) 



One proves (46) by writing out the LHS in terms of T~Lj and Q and canceling; absolute convergence 



in (48) justifies manipulations with infinite sums. □ 



7.5 The inverse of a singularly perturbed linear operator 

We will use the following immediate corollary of [AKT09, (A. 31)] 

J2 hW- ^ 2 i [ + 4 0' - !) ! ^ 6 J' ! < if 3 ^ °- 

31+32= j; ji,J2>0 



(49) 



Remark. Let us intuitively motivate the setup of the lemma 7.5 below. We start with an 



operator- valued Gevrey expansion T = Y2n>o h™^ 7 ™ which has a one-sided inverse in the form of 
the operator- valued Gevrey expansion Q = ^2 n>0 ^- n G n - We consider a third operator- valued Gevrey 
expansion T-L = ^ n> Qh n Wn. We are solving a linear equation (F + \\ a iH)u = v, where v is a known 
vector-valued Gevrey expansion and u is the Gevrey expansion to be found. We speak, in the title of 
this subsection, of a singularly perturbed linear operator because Q does not act between y s and X s 
with the same s. 



Lemma 7.5 Let a G {1, 2, ..}, C" > 1, and C > max{l, 36a a - 2 BC"e a }. 
Suppose we are given compatible with inclusions h-linear operators 



Fn,s '■ Xs — > y.s '■ 1 1 Fn^s \\x s ,y s < C" 



, n\C n 

(1- s) n ' 



0<s<l, n> 



(50) 
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and compatible with inclusions h-linear operators 

B n^C n 

forn>0, for s' > s, 9n,s',s ■ 34' -> %s \\9n,s' ,s\\ < r~t ^7T~ — 7^"- 

(s' — s) a (1 — s') n 

so that 

S~^ JO, ifN> 1; 

/ j Fj,sQm,s',s = \ id (51) 

Suppose given another family of h-linear operators compatible inclusions 

n\C n 
~Hn,s ■ X s -^y a \\H-n,s\\x a ,y a < C" ,.. - 

and elements v n = (v n ^ s ) G X\, n > 0, satisfying 

n\C n 



\v„..*\\ < 



1 - s ) n 



Then there are u n = {u n)S ) G X\ such that 



C n+a+1 (n + a + l)\ 

\ U n,s\\Xs< (1 _ s)n +a+l ' H ^° 



so that: 

If 

\h\C_ 

1-s 



< 1 (52) 



then F s = (X)„>o h n T n ,s ) ,~H S = [Yl n >oh n/ Hn,sJ ■ X s ->• y s are bounded operators, v = 
J2 n >o h n v n ,s G ys and u = J2 n >o h n u njS G X s are absolutely convergent series, and 

{F+h a U)u = v iny s . (53) 

Proof. Step 1. (Algebraic construction of the solution) 
Guided by the "formula" 

" {F+^U)- 1 = [{l + h a UF- l )F]- 1 =F- 1 {l-h a UF- 1 + h 2a {UF- l f + ...) " 

which however does not make literal sense, we take 

g{ng) n v u = " J2 ^ 0+ --- +fc " +£l+ --- +£ " +m ^o,, 1 , S o^ 1 ,, 1 ^ 1 ,, 2 , Sl ---^,s„a fcn , Sn+1 , s ^ m , Sn+1 

k , ...,k„ > 

e 1: ...,e n ,m> o 

(54) 
with sq = s and we reserve the right to choose other s^'s differently in each summand subject only to 
the condition sq < si < ... < s n+ \ < 1 (compositions will not depend on these choices). We have 

I |y&o,si,so '[■h,siyk-i_,S2,si---'Ll n ,s n yk n ,s rl+ i,s n v rn,s n+1 \\ S 
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< 



B k \C k ° C"h\C h 



(si- So ) a (l-si) fco (l-5i) £l ■"" 
B kn^lC^- 1 C"£ n \C £n B k n \C k " m\C m 

" ' (S„ - S n ^Y (1 - S„) fc -1 (1 - S n )^ ( Sn+1 - S n ) a (1 - S n+1 ) k " (1 - S n+1 ) m ' 

Choose s ra+ i so that 



•Sn+1 Sri 



m + k n + a 



(l-«n), 



then 



Then 



/ m+k n +a \g ( m+k n +a \m+kr, 
\ a / V m+k n > 



< 



(m + k n + a) a e a 



(S n+1 - S n ) a (l - S n+ i) k " +m (1 - Sn ) k «+ m+a - a a (l - Sn )k n +m+a 



RHSofd55b < } ' , w-rn-fq k \l 1 }k 1 }...£ n \k n \m\x 

4_F - a (n+l)a(l _ fl() )(n+l)a \^l-s / 



x(k n + m + a) a (fc n + m + £ n + fc n _i + 2a) a ...(/c + ••• + m + [n + l]a) a . 



But 



E 



(55) 
(56) 



k Ui\k 1 \..l n \k n \m\(k n +m,+a) a (k n +m+£ n +k n ^ 1 +2a) a ...(k +...+m+[n+l}a) a < 






< (36a° _1 ) n+1 (iV+[ra+l]a)!. 



(57) 



Indeed, on the LHS of (57) let us introduce new summation indices: m + k n =: K n > 0, 



&n + k n -i =: K n _i > 0, .... , l\ + fco =: -Ko > and use (49); let us also use that 
for any p > 0. Then 



(p+a) a 
(p+l)...(p+a) 



< a' 



a-1 



LHS of (57) < 



< (Ga^ 1 )^ 1 £ K \...K n _ x \K n \ 



K + ... + A'„ = AT; 

ic , ,.,v,, k„ > o 



(K n + a)\ {K n _ 1 +K n + 2a)\ (K + ... + K n + [n + l]o)! 
X! ' (if„-i + X„ + a)! '" (A' + ... + #„ + no)! 



< (ea - 1 )' 



E 



#<> + ... + ir„ = tf; 
K ,...,K n > 



A !...^_ 2 !^_ l!( ^ + a)! f- 1 +Kn + 2 4... {K ° + - + Kn + [n + P < 



(put ^_! := K n -\ + iir„ and use (49 1 



< (6a a 



(etc. 



E 

Ko + ... + Kn-2 + K^^N; 
K Q ,...,K n -2,K' n _ x >0 



K- I K UTC' ,o„ J^-i+gn-2 + 3a)! (JC + ... + A-;_ 2 + [n+l]a)! 



< (6a 



a—l\n+l nn 



Q n {N + [n + l]a)! < (36o a - i ) n+i (iV + [n + l]a)!, 



hence ((57). 
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Let us go back to interpreting ( |54[ ). Formally writing 

" ^2h na g(ng) n v = ^2h- 



U.I 



n>0 



J>o 



motivates the definition 



U J '■— / j 7 j ykQ,si,sa>i-l\,s 1 yki,S2,s 1 ---iLl n ,s n yk n ,s n+1 , 

na+N=J, n,N>0 k +...+e n +m=N 



n v m,s n+1 , 



where, using (57) 



Ml< E 

na+N=J, n,N>0 



( fj"\n j^n+1 e (n+l)a 



c 



N 



j(n+l)a (i _ SQ )JV+(n+l)a 



(36a a - i )" +i (iV + (n + l)a)!< 



(assuming C" > 1, C > 36a a 2 BC"e a , C > 1 and taking into account that we have < J+l < J+a + 1 
summands) 

(J + a + l)!C J+1 (J + a + l)!C J++a+1 

< ; ±r, < - 



:i - s y+ 



1 - s ) J+a+1 



Step 2. (Proof of (53)) 



Note that estimates on v n>s and u n ,s for s outside of the range prescribed by j^j < 1 will be 
used in an essential way. 

In the expression 



(T + h a n)[J2(-i) n h an J2 h ko+ ---+ k - + ^ + ---+^+ m g k0tSuso n eusl g kl ^ sl ...n in ^g kn ., Sn+1 



s n u m,s n ^.i J 



n>0 



kg, ■■■, k„ > 
£l, ...,t n , m > 



(58) 
a well-defined operator is applied to an absolutely convergent sum, once s/s j > 1, are chosen sepa- 
rately in each summand as in (56) while sq = s ; this absolute convergence justifies the manipulations 
with infinite sums below. 



Let X denote the set of all indices n, ko, .., m as in the sum (58 ); to stress dependence of choices of 



si, ..., s n+ \ on l £l we will write s\{l), ..., s n+ i(i). The entries of i will be denoted n(t), ko(i), .., m(t); 
a{i) will denote /cq + ••• + t n + m- 



Rewrite (58): 



(58) = / + //, 



where 



E (-l) n ^ an+ ' ?+ ' T(t) ^'j, s ^fco, S i(0^o^i,siW^i 1 s 2 (0, S iW 



n >0; j>0; fto, .... fcn > 
£i,...,£ n ,m > 



// 



£ 



( _ irr Ki) + M)^ 5fc0)Si(()iS0%)SiWftiiS2(()iSi(i)O .. 



">0; j > 0; fco,.-.,fcn>0 
4,...,C,m>0 
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••• ° 'T-^,i,s„(t)t'fc, l ,s n+ i((,),s n (t)' L 'm,s n+ i(t)- 



In I introduce a new index p = j + ko: 



E 



.\) n h an +P+ fe l 



fci+...+^ n +m 



X] Fj,sGk ,Si,s ° 



n>0; p>0; fei,...,fc„>0 
e u ...,£ n ,m> 



U,fco>0; i+fco=P 



°'H£i,5i^A: 1 ,s 2 (0,Si---^ n ,s„(0^A; n ,s n+1 ( t ),Sn(0 U m>«n+i(0' 

where 5i depends on n, ki, ..., k n ,£\, ...,£ n ,m,p and is denned as min{si(t)} over all is with prescribed 



n,ki, ...,k n ,£x,...,£ n ,m and ko < p (so it is a minimum over a set of p+ 1 elements). By (51) the sum 
in parentheses is zero unless p = in which case it is the inclusion y§ 1 c —t y so . Hence 



E 



(-l) n h 



nuan+ki+...+£ n +m 



%i,8^ 1 , S j,( l ), s -%, s „( l )^, s „ +1 ( 1 ), Sn ( l ) t 'm,s„ + i( t )- ( 59 ) 



7i > 0; fci, ..., k n > 
^i,...,C,m>0 



We see that most terms of I are canceled by terms of //, and only the n = term of (59) 



remains, which is just v which implies (53). □ 



8 Proof of the main result 

In this section we prove our main result Theorem |1.1[ 

The main step in the proof of theorem |1.1| is, by section [2j the solution of equation 
F S (E ,...E M -2,T) = 0, cf. @, where F s : X T S ->■ y T s , s < 1, r < r , see section |4J With re- 



spect to s, F s form a family of maps compatible with inclusions, in the sense of section |7.1| - for 
various maps below we will keep in mind and use compatibility with inclusions without always writing 



these words. The operator h of section 7.1 will be just the multiplication by h in the spaces XJ,y^, 2J, 
with \h\ = t. We will use h- linearity properties of various operators without explicitly mentioning it. 

We will take x^ ' E X[ to be the tuple (Eq, ..., Em— 2, T) of polynomials in h that solves 
F(Eq, ...,Em-2,T) up to order h bo+1 , 6q = 8, see (63), on the disc of radius po (see Notation 4.1); 



such a tuple exists by Proposition 2.1 



If we write down the expression for F s (Eq, ...,Em-2, T) as a power series in {Eq, ..., Em-2, T) G 



XT ', the series will converge provided \\hT'\ 



V'(s) 



< 1 which is definitely the case if 



r \\(Eq, ...,Em-2,T)\\xj < 1. Therefore, if we impose on r the the condition 



n Coin 1 

T \\x ( '\\x[ < 2' 



(60) 



then F s (x(°' + w) can be represented in the form (34) with some finite a independent of s. 



Since F s (Eq, ..., Em-2,T) is an analytic function of Eq, ..., Em-2,T, we can write its Frechet 
derivative by means of usual Calculus formulas. Since F s explicitly depends only on the 
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Eq, ...,Em-2,T, 

pressed) 



,T( 3 \ we can write dF( Eot ,.. t E M _ 2 ,T)(£o) ■■■,£m-2,T) as a finite sum (index s sup- 



N 



dF, 



(E , 



,T) 



(£o, ...,£m-2,T) = } y ^y(Eo, ...,Em-2,T) ■ L v (Eq, ...,Em-2,T), 



(61) 



v=l 



where (So, ...,£m-2,T) belong to the tangent space of XT at (Eq, ..., Em-2iT), and F u : XT — > yj 
are analytic on an open subset of XT compatible with inclusions XT — > XJ, and ^J — > ^J/, and 
L u : XI — > yi are constant (i.e. (Eq, ..., Em -2, T)-mdependent) linear maps XT — > yj of norm < 1, 



also compatible with inclusions. The dot (•) on the RHS of (61) is the pointwise multiplication of 



elements of yj which has norm 1 in the sense that \\y\ • j/2 1 1 3^ j < 1 1 2/1 1 1 3^J" ll^llyj- 

For concreteness, variating the LHS (13), we pick N = M + 5 and make the following choices 
where specifics of big formulas will be unimportant later on: 



Ji = 




1, 






^i=e: = o^ 


F2+3 - 


hF 2 +j - 


E' M =i (D^T^z 3 -" + hT' (2 + hT') (z 


+ 


hT) 1 


L2+3 = £3 , j 


Fm+1 = 


hFM+i 


(see below) 






Lm+i — T 


Fm+2 = 


hTm+2 


(see below) 






Lm+2 — T 


Fm+3 = 


h Tm+s = 


h 3 T'" 3h 4 T" 2 
2 {1+hT') 2 2 (1+hT') 3 






Lm+3 — T 


Fm+A = 


h J-M+4 = 


3h 3 T" 
2 (1+hT 1 ) 2 






Lm+4 — T 


Fm+5 = 


h J-M+5 = 


h 2 1 
2 1+hT' 






Lm+5 — T , 



>■ 2 ' 

0,...,M-2; 



-T; 



where 



(62) 



F M +i = (1 + hT') 2 J2 JEj(z + hT) 



3-1 



£" 2 Oh^z M ~^T^ MT > (2 + hr) {z + hT y 



3=0 



F L 



M+2 



Y,Ej(z + hT) 3 -jY. ^'^h^F 



ED 



\ 3=0 ^=1 

Notice that norms of L v are all < 1 thanks to the assumption on po that we made in Notation 



Motivated by our considerations in section 6.2, we have grouped the terms F V L U in such a way 
that terms v = 1, ..., M+2 only depend on Eq, ..., Em-2,T, T" and £q, ..., Em-2,T, T'; terms containing 
higher derivatives are put into summands for v = M+3, M +4, M+5, we notice that all summands the 
latter group also contains a factor of h 2 . We will thus treat Fm+zLm+z + Fm+^Lm+a + -^M+s-^Af+s 
as a perturbation of F\L\ + ... + Fm+2Lm+2- In turn, as F2, ■■-,Fm+2 contain a prefactor of h, we 
treat F2L2 + ... + Fm+2Lm+2 as a perturbation of F\L\. 

Modify a in such a way that representations of the form (34) also hold for F\(x^' + 

w), ...,F M +2(x {0) + w),F M+ 3(x^ + w), ...,F M+ 5(x {0) + w). 

From now on we fix the following notation: 
A is the constant such that L\ understood as a map Z T S 



y T s has an inverse of norm < A; 



B is the constant such that the restriction map ZJ 
we fix the following functions: 



XT, s' < s, is of norm < 



B 



(s-s') 2 



P(t) = te 4i ; (3 2 (t) = J2 



k>\ 



k+lAk 
k\ ' 
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We will need an integer sequence bj, j > 0, satisfying the properties: 

a) bj+i = 2bj - 7; b) b j > 3 + j; (63) 

clearly the condition a) together with bo = 8 generates such a sequence. 

We assume that r satisfies ( 60 ) , ( 67 ) , ( 71 ) , ( 73 ) , ( 79 ) , and p0| ) ; we assume that a is a finite 
number chosen as above, and we assume that C satisfies (78) and (67). 

Step 1: Constructing iterations of the Newton method as formal factorially divergent expansions. 

We will construct Gevrey expansions which will constitute iterations of the Newton method. 
For our intuition, we suggest the following correspondence (which we will not make precise) of the 
notions of section |6.1| to the objects introduced below: 



,^)^ T (0) 4-„ti) 



7/1 (J V 
w secl6H 



jH^W® and h b i- 3 ^2h n w. 



U) 



n>l 
i/secl6H 



n>l 



n>0 



Speaking rigorously, we will now construct the following objects for all j > 0: 

n\C n 



forn>o y^ = (yg)e^ r ; ||yg|| < 

forn>l w^ = (wtiheX{; ||u)g>|| < 



(1-s)"' 

1 C n n\ 
¥ (l-s) r 



We will use the notation: 



(./) 



forn>l <?(°)=0; <#> = 4°) + ... + w^, j > 1; 



7 (i) 
l5Slkr<2- 



(64) 
(65) 

(66) 



clearly #„ € -X"i and, as soon as (65) is true, <7n,s satisfy (35) with Co = 2, i.e. 

C n n! 



(!->)' 



By our choice of x*- ', we can write F(x^ ') as a convergent power series h a+1 ^2 n>0 y n (z)h n , 
where y n (z) is an /i-independent function and \\y n {z)\\y[ < d ■ c n . Put yn,s '■= hy n (z), then assuming 

C>c and r • c' < 1 (67) 



we obtain (64) for j = 0. 

Suppose y)ij are defined for j' < j and wi"| s are defined for j' < j — 1 (vacuously true if j = 0). 
Let us construct w^' and 3/W+ 1 ). 

Using lemma 7.2 with x^ ^ for xq, with g„^ for 5n, s , and with JF V and JS, for /, delivers elements 



G v , n £ 3-T, H^n.sllyj ^ /3(Co") 



C n n\ 
(l-s) n ' 
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n>l, r/ = 2,...,M + 2; 



(68) 



Gl, n ey{, WG'^Jyj < max{\\h(xW)\\yr,l3(C Q a)} ™ n > 0, v = M + 3,M + 4,M + 5, 

(69) 
such that, intuitively speaking, P v (x^ > + J2 n h n gn ) corresponds to ^(x^ ') + E n>1 ^"G Wi „, and 

A,(* (0) + En ^"5n } ) corresponds to £ n > >» n &,n- 

Let us find an infinite series of operators £ njS : yi — > Z T S so that the infinite sum En>o h n Cn,s 
will play the role of the inverse of Ylu=i ^(x) ■ L v in the precise sense specified below. An estimate 



of ||jC niS ||/ ( (i_2)n ) by 2A will be important in (77) which, in turn, lets us preserve the same constant 
C from one induction step to the other. 



With this goal in mind, we apply the lemma 7.4 with operators 



4 M+2 



(^Lemma = ^^(x^) • L v , {U n ,s)hemma. ■= ^ hG Vj n l3 ■ L v , n> 1. (70) 



v=l v=2 

and constants 



(Cl)Lemma = (M + l)r/3(C a); ^Lemma = 2A he re- 



In order to assure that, as the assumptions of lemma |7.4| require, J- has inverse of norm < 2A, we 
remember that L\ has an inverse of norm < A; it is for that reason that we assumed that r should 
satisfy 

M+2 

T . ^ \\^ x W)\\ yi < - (71) 

v=2 

We have assumed above that r is so small that 

((Ci) Lemma A) 2 e^ c ^— A < 2A, (72) 

or, more explicitly, 

((M + l)TAp{2a)fe i{M+1)TA ^ 2a) < 2A; (73) 

therefore, the lemma [7T4] yields operators 

C n n\ 



L n ,:yi^Zl s.th. \\£ n J<2A 



(l-s) n 



satisfying (46). 



For any s' < s we can compose the operator C n s with the restriction map Z T S -+ X^, and get 
£'n, s , s ' ■ yi "► X} of norm < ^ • 2A^. 



We are thus in the situation of Lemma 7.5 with a = 2, operators (cf. (|68[), (69)) 



M+2 M+2 

P^Lemma = ^ ^(a^) ' ^, (^)Lcmma = ^ G>,„ ■ L„, if 71 > 1; (74) 

y=l !/=2 

M+5 A A/+5 

(rijLemma = / •A/(aJj " i'l/ ; 1-6. (rtnJLemma = / "V,n " ^vi ('5) 
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vectors (see (64) 



(Wn)Lemma = ~Vn\ U > 0] (76) 

and constants 

M+2 

(C')Lemn» = max{ ^ \T u (x^)\ X[ , d}; 
v=l 

Af+5 

(C") Lemma = C":= Y, max{||^(x(°))b r ,/3(Coa)}; (77) 

u=M+3 
-DLcmma — ^-Dhere^; 

and we must assume 

C > max{l, 36 • 2BA ■ C" ■ e 2 }. (78) 



Then Lemma 7J5 delivers vectors u n s.th. 

C n+3 (ra + 3)! 
IKi,«||;rj < 7^ _ ^ n+3 

Put Wn = u n -3 for n > 3, u^ = wJj = 0, then u) = Xm>i ' l " ul n plays the role of the solution of 

(dF) x (w) = fc 3 J> ft i#>. 

n>l 

Define w^' = h b i~ 3 w. Because of p3| ), b) and because have assumed 

r < ^ (79) 



we see that the estimate of (65) is satisfied. 



We will now obtain vectors yfl such that, intuitively speaking, /i J+1 Xm>i h n yn plays the 
role of /(x(°) + 5^') + w&) - f(x^ + gd')) - [/'(x(°) + gCj))( w «))] (the cancellation of the second and 



third summands will be discussed later), apply lemma 7.3 b) with the following inputs: 

/Lemma — -^here; l^O/Lemma — •£ > 0Lemma — 9 i 

"Lemma — "j «j ("riJLemma — U! n , 
C^Lemma — Cohere; V^O/Lcmma — ^j ^Lemma — ■!•• 



Then Lemma 7.3 b) yields 

m«,s G 3£, n > 1 

satisfying 

n!C n n!C7 r> 

Il u n,s||yj < /32(aLemma((Co)Lemma + ^Lemma)7: r- = fa{?>Oi) 



(l-s) n y '(l-s) n 

such that h 2bj ~® X/?>i ^ u i plays the role of y" +1 \ 
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In (63), a) we have denned bj+i = lb j — 7; now take yfl = hu n and assume that 



T-[3 2 (3a)< 1. 



(80) 



,(i+i) 



Then y„ satisfy the inductive assumption (64) 



The inductive construction of step 1 is thus complete. 



Step 2. (Passing from Gevrey expansions to actual vectors) 

If condition 

tC < 1-s 



(81) 



is satisfied, then we can define the following vectors by means of absolutely convergent sums: 

r (j) = V^ h^n(J) abs -^ onv - h bj-3 ST h^,,ti) a y*. (82) 



w - - £ * n *ffi s = nv - hbj ~ 3 £ ^S e ** T ; 



n>l 



n>l 



gW) =w (0) + ... + w y-l)^^v.j-^g) G ^r. 



n>l 



y® = h h *Y,h n V$ e ^ 



(83) 

(84) 



n>0 



the equalities marked above as (abs.conv.) are justifiable by operations on absolutely convergent sums. 
Claim. If d85b and (f87|) are satisfied, then F{x® + g (i) ) = y (i) . 



The case j = is obvious by definitions. 

Suppose the Claim is true for j, let us deduce it for j + 1. 



Indeed, applying Lemma 7.3 with the same ingredients as in Step 1 but under condition 



En 



T n C r, 



< 



■* (l-s) n ' 3a 

n=l 

(which we have assumed above) corresponding to ( |42[ ), we have 

F ( x (o) + g (i) + W U)) = y (i+i) + F ( x (f>) + g (i)) + [ F '( x (o) + g 0'))]( w (i)) 

so our claim, in view of induction hypothesis, is reduced to showing that 



(85) 



rti) 



.[ F '( x (o) +g O))]( w O> 



(86) 



By Lemma 7.5, w'-" solves the equation 

(F + /i 2 H)w^ = -y® 



where (cf. d74b, (75) 



f = e^(* (0) ) • L » +X> n X>*,», s • L - 



v=\ 



n>\ v=2 
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H — 2_j h n /_^ G uns ■ L v 



n>0 



v=5 



But the definition of G u and G' v by means of lemma 7.2 shows that under condition 

1 



^ r n C n 
^ (l-s) n < 

n>l v ' 



2a 



(87) 



(which we have assumed) corresponding to ( 36 ) , we have 

4 7 

and so by ( |ol~] ) we conclude the proof of the claim. 

Step 3. (Proving the convergence of the Newton's method.) 



By the above construction, 



n>0 V ' 



, 1 + 

2J V 3a 



w 



U)\ 



XI 



< 



-y- 

21 ^ (1 



T n C r 



<(l-8)' 



< 



Thus x^ + g^ has a limit a^ ) + 



g 



in XI and F(x(°) + g°°) 



3a2^ 



0. 



Let us review the choices of various constants that we have made. In the beginning of this 
section we have chosen and fixed a large enough depending on the initial data (namely, on Qi(z,h)). 



In Step 1 we have assumed that r is small enough to satisfy (60), (67), (71), (73), (79), and 



(which depend only the initial data); then we have chosen and fixed C large enough to satisfy 



and (67); a choice of C valid for one value of r is also valid for smaller values of r. In Step 2, in 



80 



we choose and fix s; a choice of s valid for one value of r will also be valid for smaller values of r. 



Then we shrink r if necessary to satisfy (85) and (87). We conclude that there exist s* and r* such 



that for any s < s* and r < r*, the equation F(x) = has a solution x G XJ . 



Finally, in view of sectional we have shown that Eq, ...,Em-2, T, and hence also y(x,h) have 
the Gevrey growth condition. Using (ml), one shows by a standard argument that also ip(x, h) satisfies 
a Gevrey growth condition. This concludes the proof of theorem |1.1| □ 
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